Layer-number determination in graphene by out-of-plane phonons 
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We present and discuss a double-resonant Raman mode in few-layer graphene, which has not 
been interpreted before and is able to probe the number of graphene layers. This so-called A'^ mode 
on the low-frequency side of the G mode results from a double-resonant Stokes/anti-Stokes process 
combining an optical (LO) and an out-of-plane (ZO') phonon. Simulations of the double-resonant 
Raman spectra in bilayer graphene show very good agreement with the experiments. 
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Raman spectroscopy belongs to the most widely used 
methods in graphene research. Raman spectroscopy is 
used for characterizing graphene regarding d efect a^ ^, 
doping^, strairpnS, crystallographic oricntatiorpSHll^ or 
interaction with the substrat^^^. In view of fundamen- 
tal physical properties of graphene, Raman spectroscopy 
gives information on electron-phonon coupling and scat- 
tering rates, optical excitations in graphene, thermal and 
mechanical propertiea^'^ ^^. Probably the most popular 
application of Raman scattering in graphene is the dis- 
tinction of single-layer graphene from few-layer graphene 
and graphite via the lineshape of the double-resonant 2D 
models. On the other hand, few-layer graphene has re- 
cently come into focus, a s gate d bi- and trilayer graphene 
offer a tunable band gap F^^^ I and bilayer graphene has 
been demonstrated to give much higher on-off ratios 
in a field-effect transistor than single-layer graphene^^. 
Therefore, it is important to establish a reliable method 
for the determination of the layer number in few-layer 
graphene and to identify spectroscopic signatures of the 
layer-layer interaction. So far, typically the evolution of 
the 2I?-mode lineshape or the absolute Raman intensity 
of the G mode is used in combination with optical con- 
trast measurements. However, the lineshape of the 2D 
mode depends strongly on the excitation wavelengtniS^ 
and the G-modc amplitude depends not only on the scat- 
tering volum^2012ij^ ^^^ g^jgQ Qjj ^jjg substrate and opti- 
cal interference effects^^l. Recently, the rigid-layer shear 
mode, which is the other Raman-active E2g phonon mode 
in graphite, was shown to have a strong frequency depen- 
dence on the number of layers in few-layer graphene^^l. 
The frequency of this mode, however, is below 44 cm^^. 
Measurement of this low-frequency mode is therefore dif- 
ficult and requires non-standard equipment. 

Here we present and interpret a newly discovered Ra- 
man mode on the low- frequency side of the G mode, 
which can be used to determine the number of layers 
in few-layer graphene. This so-called N mode is based 
on a double-resonant intravalley scattering process com- 
bining the longitudinal optical (LO) and the rigid-layer 
compression mode (ZO'). The peak position as well as 
the lineshape of this peak allow an assignment of the Ra- 
man spectra to the number of graphene layers for up to 
approximately eight layers. In addition, we simulate the 
double-resonant Raman spectra in the iV-mode region for 
various excitation energies in bilayer graphene. 
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FIG. 1. (a) Ratio of Color Difference (RCD) measurement on 
a few-layer graphene sample, (b) Few-layer graphene and the 
scanning path on the sample. Lengths along the RCD scan 
are not drawn to scale. 



Graphene samples were prepared by mechanical cleav- 
age under cleanroom conditions from natural graphite 
flakes and transferred onto a silicon substrate with an ox- 
ide thickness of 80 nm. The samples were analyzed with 
an optical microscope (Olympus BX51M with an 100 x 
objective) . We determined the number of graphene layers 
by optical contrast, using the Ratio of Color Difference 
(RCD) method^''. The RCD values were calculated using 
the formalism from Ref.[25] 
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where Xq, Yq, and Zq denote the tristimulus color com- 
ponents of the Si/Si02 substrate and X„, !"„, and Z„ the 
color components of n-layer graphene. Since the RCD is 
independent of the light source^Sl, the RCD values can be 
calculated directly from the RGB color values of the opti- 
cal image. In Fig.flTa) an exemplary result of a RCD scan 
is shown. Here, the RCD measurement along the high- 
lighted path revelled graphene thicknesses ranging from 
71=2 to n=6 layers. This result corresponds to the opti- 
cal contrast from the image, which is shown in Fig.mb). 
Graphene samples with layer thicknesses up to 11 lay- 
ers have been prepared and were characterized by this 
method. 

We performed confocal /x-Raman measurements under 
ambient conditions using a LabRAM HR800 spectrom- 
eter. Laser excitation wavelengths of 532 nm (2.33 eV) 
and 633 nm (1.96 eV) were chosen. Raman spectra were 
recorded in back-scattering geometry with a spectral res- 
olution better than 1 cm"-*^. The laser was focussed with 




1460 1480 1500 1520 1540 1560 1580 
Raman Shift fcm^''^) 



FIG. 2. Raman spectra of n-layer graphene on the low- 
frequency side of the G mode at 633 nm laser wavelength. 
Spectra are normalized to the same G-mode amplitude and 
vertically offset for clarity. The dashed lines serve as a guide 
to the eye. 



an 100 X objective and had a spot size <500nm. All spec- 
tra were calibrated by standard atomic emission lines of 
neon (Ne). 

The bandstructure and phonon dispersion of bilayer 
graphene were calculated using the Siesta DFT code 
in local-density approximation^''. The calculations were 
performed according to Ref.l^H We used the exper- 
imental geometrical values of graphite, i.e., a lattice 
constant of a = 2.46 A and an interlayer distance of 
c/2 = 3.35 A12S1. The T-point frequency of the £2^^ mode 
was scaled by a factor of 0.96 to the experimental value 
of 1584 cm^^. We rescaled the calculated phonon disper- 
sion by the same factor; the resulting phonon dispersion 
shows very good agreement with experimental datct^. 

Fig. [2] shows the Raman spectra of n-layer graphene 
for layer thicknesses ranging from monolayer to 11-layer 
graphene at 633 nm laser wavelength. For 71 > 2 we ob- 
serve a layer-dependent peak on the low-frequency side 
of the G mode. This mode is approximately 100 times 
weaker than the G mode. It is clearly absent in mono- 
layer graphene. Furthermore, additional peaks appear 
for more layers. We label these Raman modes in the 
order of their appearance as iVi, N2, and N^ and refer 
to them as N mode. The layer-dependent shift of their 
peak position is shown in Fig. [3] The frequencies decrease 
and tend towards a lower limit as the layer thickness is 
increased. 

The absence of the N mode in single layer graphene 
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FIG. 3. Peak positions of the A^ mode as a function of the 
number of graphene layers at 633 nm laser wavelength. 



indicates that it may originate from interlayer vibra- 
tions. We assign the N mode to a double-resonant in- 
travalley scattering close to the K point combining LO 
(longitudinal optical) and ZO' (rigid-layer compression) 
phonons, in which the LO phonon is Stokes- and the ZO' 
phonon anti-Stokes scattered. An illustration of the dou- 
ble resonance is shown in Fig.|4] The dashed horizontal 
line corresponds to the phonon frequency of the defect- 
scattered LO phonon, i.e., the D' mode. In (a) the elec- 
tron is first scattered by an LO phonon and afterwards 
the Stokes/anti-Stokes scattering with a ZO' phonon fol- 
lows. The reversed order in the scattering process is 
shown in Fig.WTb). The double-resonant scattering in 
a two-dimensional illustration is shown in Fig.HTc). The 
resonantly enhanced phonon wave vector along T-M con- 
nects two electronic states on the K-M high-symmetry 
line. An explanation of this scattering process is given 
below. 

Our assumption is supported by the correspondence 
between the N mode and the double-resonant LO-I-ZO' 
peak (~1740cm~^), resulting from an intraval ley d ouble 
resonance combining an LO and ZO' phonorpSEfll ^g 
label this peak in the following as LOZO'"*". The Raman 
spectra of the LOZO'"*" peak for layer thicknesses from 
monolayer to 11-layer graphene are shown in Fig. [5] All 
peaks of the N mode and the LOZO'"*" peak have ap- 
proximately the same distance to the D' mode, which 
can be resolved at ~1616cm~^ for 633 nm laser wave- 
length. Due to this symmetry, both the N mode and 
LOZO'+ peak must differ from the D' mode in the same 
process, namely by the scattering with a ZO' phonon. In 
the case of the N mode, the ZO' phonon is anti-Stokes 
scattered, whereas the ZO' phonon is Stokes-scattered for 
the LOZO'"*" peak. This combination of Stokes and anti- 
Stokes scattered phonons in a double-resonant process 
was never reported before. 

Since we assign the N mode to a double-resonant Ra- 
man process, we would expect a laser energy depen- 
dent shift of the peak position, as this is a fingerprint 
of double-resonant Raman scattering. Fig.|6] shows the 
spectra of bilayer graphene for 532 nm and 633 nm exci- 
tation wavelength. The LOZO'+ peak blueshifts with 
increasing laser wavelength, whereas a shift of the N 





FIG. 4. Schematic illustration of the double-resonant in- 
travalley scattering of a LO phonon (blue, dotted) and a ZO' 
phonon (green, dashed). Energies are not drawn to scale. In 

(a) the electron is first scattered by the LO phonon, then 
the Stokes/anti-Stokes scattering with a ZO' phonon follows. 

(b) shows the reversed order, (c) Contour plot of graphene's 
bandstructure around a K point. The K-M direction is high- 
lighted in white. The orange arrow denotes the phonon wave 
vector q that is enhanced by the double resonance. 
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FIG. 5. Raman spectra of the double-resonant LOZO'^ peak 
in n-layer graphene at 633 nm laser excitation wavelength. 
Spectra are normalized and vertically offset for clarity. The 
dashed lines serve as a guide to the eye. 



mode cannot be observed or is on the order of our spec- 
tral resolution. This behavior can be understood from 
the dispersion of the LO and ZO' phonon branch shown 
in Fig.lWa). When a double-resonant process combines 
Stokes and anti-Stokes scattered phonons, the difference 
of both phonon frequencies determines the final peak po- 
sition. The LO and ZO' phonon branch exhibit nearly the 
same slope in the relevant range. Hence, the difference 
of both phonon branches is nearly constant. Therefore, 
a change of the phonon wavevector does not result in a 
shift of the phonon frequency and the N mode shows no 
or little dispersion in the range between 1.9 eV and 2.3 eV 
laser energy. In fact, the shift of the TV mode between 
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FIG. 6. Comparison of the N mode and LOZO^ peak in bi- 
layer graphene for 532 nm (blue, lower spectrum) and 633 nm 
(black, upper spectrum) laser excitation wavelength. 



633 nm and 532 nm excitation wavelength, derived from 
the phonon dispersion in Fig.^a), is less than lcm~^. 
This result fits our experimental observations very well. 
For the LOZO'''" peak, both phonon branches must be 
added. The resulting phonon branch has a positive slope; 
therefore the peak position should increase for higher ex- 
citation energies. We estimate from Fig.JTJJa) a blueshift 
of ~7cm~^, which is close to the experimentally obtained 
shift of -1-5 cm^^. 

The shift of the Ni peaks [i — 1,2,3) as a function 
of the number of graphene layers can be explained with 
the evolution of the ZO' phonon spectra in few-layer 
graphene. In n-layer graphene there ex ist n- 1 vibra- 
tions with rigid-layer compression patterrplESl. There- 
fore, for an increasing number of graphene layers, the 
LO phonon can scatter with an increasing number of ZO' 
phonons. The ZO' vibrations exhibit a layer-dependent 
shift towards an upper limit, i.e., the frequency in bulk 
graphitc'^^ ^^. This explains the downshift of the N mode 
and the upshift of the LOZO'+ peak as a function of the 
number of graphene layers, as well as the appearance of 
additional peaks for increasing number of layers. Recent 
work from Lui et al. show a similar layer-dependence of 
the LOZO'+ mode on the high-frequency side of the G 
mode for up to six layers^, in agreement with the spectra 
shown in Fig.[5J 

To support our interpretation, we simulated the 
double-resonant Raman spectra using the equatioriS^ 
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where El is the energy of the incoming photon and M 
are the matrix elements, which are assumed to be con- 
stant. However, the strong angular dependence of the 
optical matrix elements was taken into account by set- 
ting the integration path as shown in Fig.[4£c), in agree- 
ment with results for the D' mode in Re™. Thus, the 
optical transitions are calculated along K-M, whereas 
the phonons predominantly stem from the T-M direc- 
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FIG. 7. (a) LO phonon branch (dashed), LO-ZO phonon 
branch (lower curve), and the LO+ZO phonon branch (upper 
curve) along the F-M direction in bilayer graphene. The M 
point corresponds to 2 7r/(\/3a) with a = 2.46 A. The ver- 
tical dotted lines denote the phonon wavevector at 1.96 eV 
and 2.33 eV laser excitation energy, respectively, (b) Sim- 
ulated double-resonant Raman spectra in bilayer graphene 
(blue curve) in the A'' mode region compared to experimen- 
tal spectra (black curve), (c) Laser energy dependent peak 
positions of the N mode calculated for bilayer graphene. 



Results of our calculations for bilayer graphene are 
shown in Fig.^b) and (c). The simulated spectra in 
Fig.IWb) fit our experimental data very well. The laser- 
energy dependent peak position of the N mode is shown 
in Fig.JTlc). The N mode follows the dispersion of the 
LO-ZO phonon branch. The laser-dependent peak shift 
is in the visible range much less than that of the LOZO"*" 
peak, in agreement with the experiments. At higher ex- 
citation energies above ~2.5 eV, we observe a splitting of 
the TV mode due to distinct contributions from the two 
IT bands in bilayer graphene. 

In summary, we have presented and interpreted a layer- 
number dependent Raman mode on the low-frequency 
side of the G mode in few-layer graphene. This so-called 
N mode is a combination mode of a Stokes-scattered LO 
phonon and an anti-Stokes scattered ZO' phonon. The 
investigation of the peak positions enables determination 
of the number n of graphene layers up to n = 8. 

The simulation of the double-resonant Raman spec- 
tra agrees very well with the experimental results. The 
N mode shows in the visible range only little dispersion 
with laser wavelength. Furthermore, the N mode does 
not overlap with other overtones or combinational modes, 
in contrast to the LO+ZO' peak. Depending on the exci- 
tation wavelength, the N mode may also be indicative of 
the stacking order in few-layer graphene. Furthermore, 
the study of ZO' phonons can give information about the 
strenght of layer-layer interaction in few-layer graphene. 

Since the occurence of the ZO' vibration is not re- 
stricted to graphene, this approach of determining the 
number of layers might be transferable to other layered 
materials. 



tion. The energy differences between the intermediate 
electronic states a, b, c and the initial state i are labeled 
as Exi- The broadening factor 7 was set to 40meV^. 
Our calculations also include the reversed order, where 
the ZO' phonon is scattered first and the LO phonon 
secondly, and scattering by both electrons and holes. 
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